Prospects of high temperature ferromagnetism in (Ga,Mn)As semiconductors 
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We report on a comprehensive combined experimental and theoretical study of Curie temperature 
trends in (Ga,Mn)As ferromagnetic semiconductors. Broad agreement between theoretical expec- 
tations and measured data allows us to conclude that T c in high-quality metallic samples increases 
linearly with the number of uncompensated local moments on Mn&i acceptors, with no sign of 
saturation. Room temperature ferromagnetism is expected for a 10% concentration of these local 
moments. Our magnetotransport and magnetization data are consistnent with the picture in which 
Mn impurities incorporated during growth at interstitial Mni positions act as double-donors and 
compensate neighboring MnGa local moments because of strong near-neighbor MnGa-Mm antiferro- 
magnetic coupling. These defects can be efficiently removed by post-growth annealing. Our analysis 
suggests that there is no fundamental obstacle to substitutional MnGa doping in high-quality ma- 
terials beyond our current maximum level of 6.2%, although this achievement will require further 
advances in growth condition control. Modest charge compensation does not limit the maximum 
Curie temperature possible in ferromagnetic semiconductors based on (Ga,Mn)As. 

PACS numbers: 75.50.Pp,75.30.Gw,73.61.Ey 



I. INTRODUCTION 

After some frustration in the community caused by 
the difficulties encountered in overcoming the apparent 
Curie temperature limit in (Ga,Mn)As of T c = 110 k. 1 2 1 
the record transition temperature has been steadily in- 
creasing over the last two years A^LS The maximum 
T c = 173 K reported 8 - to date is likely another short- 
lived record in bulk (Ga,Mn)As ferromagnets. It is now 
established that the success has been made possible by 
the technological progress in controlling crystallographic 
quality of the materials, namely, in reducing the number 
of unintentional charge and moment compensating de- 
fects through optimized growth and post-growth anneal- 
ing procedures' 7 ^ Experiments also suggest that the 
general picture of ferromagnetism that applies to these 
metallic (Ga,Mn)As systems is the one in which magnetic 
coupling between local Mn moments is mediated by de- 
localized holes in the (Ga,Mn)As valence band. The fact 
that the mechanism does not imply a fundamental T c 
limit below room temperature motivates a detailed anal- 
ysis of our understanding of the T c trends in currently 
available high quality metallic materials with Mn doping 
ranging from approximately 2% to 9%. 

Curie temperatures in metallic (Ga,Mn)As 
have been studied theoretically starting from 
semi-phenomenologicati2i 1 ii2ii&i 4 i and microscopic 
modelo 15 ' 16 i 17 i 18 i 19 ' 20 of the electronic structure. The 
former approach asserts a localized character of the 
five MnGa d-orbitals forming a moment S — 5/2 and 



describes hole states in the valence band using the 
Kohn-Luttinger parameterization for GaAs^i and a 
single constant J p d which characterizes the exchange 
interaction between MnQ a and hole spins. The exchange 
interaction follows from hybridization between Mn 
d orbitals and valence band p orbitals. The semi- 
phenomenological Hamiltonian implicitly assumes that 
a canonical transformation has been performed which 
eliminated the hybridization. 20 In this approach the 
hybridization is implicitly assumed to be weak in several 
different ways, and the canonical transformation ignored 
in representing observables. Although this approach 
is consistent, it should be realized that the localized 
d-orbitals in the phcnomcnological Hamiltonian are in 
reality hybridized with the valence band. 

The advantage of the semi-phenomenological approach 
is that it uses the experimental valu o 22 i 23 for J p d — 
54 ± 9 meV nm 3 , i.e., it correctly captures the strength 
of the magnetic interaction that has been established 
to play the key role in ferromagnetism in (Ga,Mn)As. 
The model also accounts for strong spin-orbit interac- 
tion present in the host valence band which splits the 
three p-bands into a heavy-hole, light-hole, and a split- 
off band with different dispersions. The spin-orbit cou- 
pling is not only responsible for a number of distinct 
magnetio 24 i 25 i 26 i 27 and magneto-transport 28 i 29 i 30 i 31 prop- 
erties of (Ga,Mn)As ferromagnets but the resulting com- 
plexity of the valence band was showni 3 ^ to play also an 
important role in suppressing magnetization fluctuation 
effects and, therefore, stabilizing the ferromagnetic state 
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itself. Describing the potentially complex behavior of 
MnQa in GaAs by a single parameter may oversimplify 
the problem. The calculations omit, for example, the 
contribution of direct antiferromagnetic superexchangc 
to the coupling of near-neighbor Mn pairs, and the whole 
model inevitably breaks down if valence fluctuations of 
MnQa d-electrons become strong. 

Microscopic theories, whether based on the tight- 
binding-approximation (TBA) parameterization of ener- 
gies and overlaps of valence orbitals in the lattic e) 1 ^ 21 ! 1 
or on the ab initio local-density-approximation (LDA) 
schemes ji^*iSiiLi& make no assumption upon the char- 
acter of MnQa impurities in GaAs and their magnetic 
coupling. They are therefore useful for studying material 
trends in T c as a function of Mn doping or density and po- 
sition in the lattice of other intentional or unintentional 
impurities present in real systems^ Because spin-orbit 
interactions add to the numerical complexity of calcula- 
tions that are already challenging, they have normally 
been neglected. Another shortcoming of the ab initio ap- 
proaches is the incomplete elimination of self-interaction 
effects which leads to smaller relative displacement of the 
Mn d-lcvcls and the top of the valence band. This results 
in an overestimated strength of the p-d exchange as com- 
pared to experiment. 

Within the mean-field approximation, which consid- 
ers thermodynamics of an isolated Mn moment in an ef- 
fective field and neglects correlated Mn-Mn fluctuations, 
microscopic calculations 1 - typically yield larger T c 's than 
the semi-phenomenological modelsi 2 ^ that use the ex- 
perimental J p d value. Stronger p-d exchange in the mi- 
croscopic theories leads, however, also to a larger suppres- 
sion of the Curie temperature due to fluctuation effects, 
especially so in highly-doped systems H (A closer agree- 
ment in the character of the T c versus Mn-doping curves, 
calculated within the two formalisms, is obtained when 
the deficiencies of density-functional theories are partly 
eliminated by introducing a correlation energy constant 
within the LDA+U schemes^) Despite the above weak- 
nesses of semi-phenomenological and microscopic cal- 
culations, an overall, qualitatively consistent picture is 
clearly emerging from these complementary theoretical 
approaches that, as we discuss below, provides a useful 
framework for analyzing measured T c 's. 

In experimental Curie temperature studies it is cru- 
cial to decouple intrinsic properties of (Ga,Mn)As fer- 
romagnets from extrinsic effects due to the presence of 
unintentional impurities. Arsenic antisites (Asca) and 
interstitial manganese (Mni) represent two major sources 
of charge compensation in (Ga,Mn)As grown by low- 
temperature molecular beam epitaxy (LT-MBE), both 
acting as double-donors^^. A Mni cation when at- 
tracted to a MnQa anion compensates also the MnQ a lo- 
cal moment as the two species are expected to couple 
antiferromagneticallji^SiSI due to superexchange over 
the whole range from strong to weak charge compensa- 
tion. 

The Asca antisites are stable^ up to ~ 450° C which is 



well above the transition temperature from a uniform di- 
luted magnetic semiconductor to a multiphase structure 
with metallic MnAs and other precipitates. Therefore, 
the number of Asca defects has to be minimized already 
during the LT-MBE growth by precisely controlling the 
stoichiometry of deposited epilayers^ The Mni impu- 
rity concentration can be significant in as-grown struc- 
tures. These defects are, however, much more mobile 
than the As antisites. During annealingi^^S at temper- 
atures close to the MBE growth temperature ~ 200°C 
they out-diffuse and are passivated at the epilayer sur- 
face. In this paper we have collected data for a set of 
samples that show very weak charge and moment com- 
pensation after annealing, i.e. a negligible number of 
Asca, which allows us to determine experimentally T c 
trends related to intrinsic properties of (Ga,Mn)As fer- 
romagnets. 

The paper is organized as follows: In the theory Sec- 
tion [H] we start with the semi-phenomenological mean- 
field approximation fSection lll A|) to set up a scale of ex- 
pected Curie temperatures in the material, assuming ho- 
mogeneous distribution of Mno a ions (the virtual-crystal 
approximation). We then discuss various physically dis- 
tinct effects that are not captured by this picture. In Sec- 
tion lll Bl wc evaluate the Stoner enhancement of the Curie 
temperature due to hole-hole exchange interaction. Sup- 
pression of T c due to antiferromagnetic superexchange 
contribution to the near-neighbor MnQ a -MnGa coupling 
in highly compensated samples^ is illustrated in Sec- 
tion III CI In this section we discuss also effects on T c 
arising from the discreteness of random MnQ a positions 
in the lattice that becomes important in the opposite 
regime, i.e., in systems with low charge compensation or 
co-doped with additional non-magnetic acceptors. Ef- 
fects beyond the mean-field approximation, namely the 
disappearance of the ferromagnetic long-range order due 
to collective Mnc a moments fluctuations are discussed in 
Section I'll Dl Since the Mnc a , Mni, and hole densities 
represent key parameters in the discussion of measured 
Curie temperatures, we present in Section III El theoret- 
ical predictions for equilibrium partial concentrations of 
substitutional Mnca and interstitial Mni impurities in 
as-grown samples, and in Section III Fl we estimate the 
accuracy of the Hall measurement of hole density in the 
polarized (Ga,Mn)As valence bands. 

Measured T c and hole densities are presented in Sec- 
tion IIII Al for a set of samples with different nominal Mn- 
doping, before and after annealing. Motivated by the 
above theoretical analysis we determine in Section fill Bl 
the partial density of MnQ a and Mni, and the effective 
density of uncompensated MnQ a local moments in our 
samples. The interpretation is based on total Mn-doping 
values, obtained from secondary ion mass spectroscopy 
(SIMS), and Hall measurements of the hole densities be- 
fore and after annealing. Consistency of these results is 
checked by comparisons with independent magnetization 
measurements. In Section IIII CI we present experimen- 
tal T c dependencies on uncompensated Mnc a moment 
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and hole densities and compare the data with theory 
predictions. Technological issues related to the growth 
of (Ga,Mn)As epilayers with large Mn concentrations 
are discussed in Section IIVI Our perspective on high- 
temperature ferromagnetism in (Ga,Mn)As semiconduc- 
tors is summarized in Section Ivl 



II. THEORY 



Qualitative implications of this T c -equation can be 
understood within a model itinerant hole system with 
a single spin-split band and an effective mass to*. The 
kinetic energy contribution, e^, to the total energy of the 
band holes gives a susceptibility: 



d 2 e k 

Xf ' k ~ dh 2 MF ~ 4ir 2 h 2 



(5) 



A. Mean-field virtual crystal approximation 

The description of ordered states in (Ga,Mn)As is 
greatly simplified by the virtual-crystal approximation 
in which the random Mnc a distribution is replaced by 
a continuum with the same average local moment den- 
sity and the role of other defects is neglected, apart from 
the potential hole or moment compensationiifliii^ 2 ^ Mi- 
croscopic TBA calculations showed^ very little effect of 
positional disorder on the strength of magnetic couplings 
in (Ga,Mn)As epilayers with metallic conductivities of 
interest here, which partly justifies the virtual-crystal 
approach. Other detailed theoretical studies, corrobo- 
rated by experimental data below, confirm the absence of 
any significant magnetic frustration in this ferromagnetic 
semiconductor associated with the random positions of 
MiiQa moments in the lattice^^ 

In the mean-field approximationiSiii, each local MnQa 
moment is described by a Hamiltonian Si -Hmf where Si 
is the MiiQa local spin operator, Hmf = Jpd(s), and (s) 
is the mean spin density of the valence band holes. Hmf 
is an effective field seen by the local moments due to spin- 
polarization of the band holes, analogous to the nuclear 
Knight shift. Similarly Hmf = J p d,NMn(S) is an effec- 
tive magnetic field experienced by the valence band holes, 
where (S) is the mean spin polarization of the Mno a lo- 
cal moments, and -/Vmii = 4x/af c is the Mno a density 
in Gai-^Mn^As with a lattice constant a\ c . The depen- 
dence of (S) on temperature and field Hmf is giveni 3 - by 
the Brillouin function: 



(S) = 



H 



MF 



SB s (S\H MF \/k B T) 



(1) 



I -Hmf | 

The Curie temperature is found by linearizing Hmf an d 
B s around (S) = 0: 

Hmf ~ Jpd^Mn(S)xf 

u ~ 5+1 S\H M f\ /„s 

Bs ~ — q 1 r ■ 

3 k B lc 

Here Xf is the itinerant hole spin susceptibility given by 
d(s) d 2 er 



Xf 



dh 



MF 



dh 2 

MF 



(3) 



and ex is the total energy per volume of the holes. 
Eqs. and © give 



, rp N Mn S(S+l) t2 
k B T c = J pd Xf 



(4) 



where kp is the Fermi wavevector. Within this approx- 
imation T c is proportional to the MnQ a density, to the 
hole Fermi wavevector, i.e. to p 1 / 3 where p is the hole 
density, and to the hole effective mass m*. 

To obtain quantitative predictions for the critical tem- 
perature, it is necessary to evaluate the itinerant hole 
susceptibility using a realistic band Hamiltonian, H = 
Hkl + s- fiMF, where Hkl the six-band Kohn-Luttinger 
model of the GaAs host bandSi and s is the hole spin 



operator 



11.24,25 



The results, represented by the solid 



black line in Fig. ^ are consistent with the qualita- 
tive analysis based on the parabolic band model, i.e, 
T c follows roughly the ~ xp 1 ^ 3 dependence. Based on 
these calculations, room temperature ferromagnetism in 
(Ga,Mn)As is expected for 10% Mnc a doping in weakly 
compensated samples. 



B. Stoner enhancement of T c 

In highly doped (Ga,Mn)As epilayers the hole-hole cor- 
relation effects are weak and can be neglected. The ex- 
change total energy e x adds a contribution to the hole 
spin susceptibility: 



d 2 e x 

un MF 



(6) 



which for a single parabolic spin-split band reads, 



XU 



e 2 {m*) 2 
47r 3 eft 4 



(7) 



where e is the dielectric constant of the host semiconduc- 
tor. Eq. JJJ suggests a hole-density independent Stoner 
enhancement of T c proportional to MnQa concentration 
and (to*) 2 . 

As in the non-interacting hole case discussed above, 
the detailed valence-band structure has to be accounted 
for to make quantitative estimates of the Stoner T c en- 
hancement. The red line in Fig. ^ shows the Stoner T c 
enhancement calculated numerically from Eqs. ©. As 
expected, T c stays roughly proportional to xp 1 ^ 3 even if 
hole-hole exchange interactions are included, and the en- 
hancement of the Curie temperature due to interactions 
is of order - 10 - 20%. 
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FIG. 1: Ferromagnetic transition temperatures of (Ga,Mn)As 
calculated within the effective Hamiltonian and virtual crystal 
approximation: mean- field (thick black line) , Stoner enhance- 
ment of T c (thin red line), spin- wave suppression of T c (blue 
symbols) . 



C. Discreteness of random Mnca positions, 
superexchange 

So far, the mean-field analysis of T c has neglected 
discreteness in the random MnQ a positions in the lat- 
tice and other magnetic coupling mechanisms, beside 
the kinetic-exchange, in particular the near neighbor su- 
perexchange. The former point can be expected to affect 
T c at large hole densities, i.e., when the hole Fermi wave- 
length approaches inter-atomic distances, fn the opposite 
limit of strongly compensated systems, where the over- 
all magnitude of the hole-mediated exchange is weaker, 
antiferromagnetic superexchange can dominate the near- 
neighbor Mnca-Mnca coupling^ leading to a reduced 
Curie temperature. 17,1 - This type of magnetic interac- 
tion was ignored in the previous section. We emphasize 
that the phenomenological model cannot be applied con- 
sistently when nearest neighbor interactions dominate, 
since it implicitly assumes that all length scales are longer 
than a lattice constant. We also note that net antiferro- 
magnetic coupling of near-neighbor MnG a -MnGa pairs is 
expected only in systems with large charge compensa- 
tions. In weakly compensated (Ga,Mn)As the ferromag- 
netic contribution takes over4L2& 

Beside the above effects of random Mn distribution, 
Mn positional disorder can directly modify p-d interac- 
tion when the coherence of Bloch states becomes sig- 
nificantly disturbed. Microscopic theories, such as the 
TB A/CPA calculations^ presented in this section or ap- 
proaches based on ab-initio LDA band structure^ 7 -, cap- 
ture all these effects on an equal footing and can be used 
to estimate trends in mean-field T c beyond the virtual 



crystal approximation. The theories do not assert any 
specific magnetic coupling mechanism from the outset. 
Instead, these follow from the microscopic description of 
the electronic structure of the doped crystal. 

Within the CPA, disorder effects appear in the finite 
spectral width of hole quasiparticle states. Since realiza- 
tions with near-neighbor MnQ a ions are included within 
the disorder-averaged TB A/CPA with the proper statis- 
tical probability, short-range local moment interactions 
(such as superexchange) contribute to the final magnetic 
state. 

The parameterization of our TBA Hamiltonian was 
chosen to provide the correct band gap for a pure GaAs 
crystal^ 7 , and the appropriate exchange splitting of the 
Mn d-states. Local changes of the crystal potential at 
MiiQa, represented by shifted atomic levels, are estimated 
using Ref.^3- Long-range tails of the impurity potentials, 
which become less important with increasing level of dop- 
ing, are neglected. (Note, that the Thomas-Fermi screen- 
ing length is only 3-5A for typical carrier densities^, i.e., 
comparable to the lattice constant.) Also lattice relax- 
ation effects are neglected within the CPA. 

In our TBA/CPA calculations, hole density is varied 
independently of Mnc a doping by adding non-magnetic 
donors (Si or Se) or acceptors (C or Be). The result- 
ing valence-band splitting is almost independent of the 
density of these non-magnetic impurities at fixed Nmu, 
which indicates that quasiparticle broadening due to po- 
sitional disorder has only a weak effect on the strength of 
the kinetic-exchange coupling. We intentionally did not 
use Mni donors in these calculations to avoid mixing of 
the (arguably) most important effect of this defect which 
is moment compensation. This is discussed separately 
below in Section III El 

The TBA/CPA Curie temperatures are obtained us- 
ing the compatibility of the model with the Weiss mean- 
field theory. The strength of the Mnoa-Mnca coupling is 
characterized by the energy cost of flipping one MnQ a 
moment, which can be calculated for a given chemi- 
cal composition. 50 This effective field H e ff corresponds 
to Hmf in the semi-phenomenological kinetic-exchange 
model used in the previous section, i.e., 



k B T c = —H e ff . (8) 

In Fig. |U we plot the mean-field TBA/CPA transition 
temperatures as a function of hole densities for several 
MiiQa concentrations. Since the typical T c 's here are sim- 
ilar to those in Fig. ^ we can identify, based on the com- 
parison between the two figures, the main physical origins 
of the deviations from the T c ~ xp 1 ^ 3 trend. Black dots 
in the left panel of Fig. which correspond to a rela- 
tively low local MnQa moment concentration (x = 2%) 
and hole densities ranging up to p — ANmu show the 
expected supression of T c at large p. The effect of su- 
perexchange in the opposite limit is clearly seen when 
inspecting, e.g., the x = 10% data for p < 1 nm -3 . The 
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mean- field TBA/CPA Curie temperature is largely sup- 
pressed here or even negative, meaning that the ferro- 
magnetic state becomes unstable due to the short-range 
antifcrromagnetic coupling. Also, the inhomogeneity of 
the carrier distribution in the disordered mixed crystal 
may contribute to the steep decrease of T c with increas- 
ing compensation. Although the Curie temperatures in 
the left panel of Fig. [21 appear to depart strongly for the 
T c ~ xp 1 / 3 dependence, the linearity with x is almost 
fully recovered when T c is plotted as a function of the 
number of holes per MnQ a , vl^Mn (see right panel of 
Fig.|2J). Note that for compensations (1 — p/NMn) reach- 
ing 100% this property of the superexchange coupling is 
reminiscent of the behavior of (II,Mn)VI diluted mag- 
netic semiconductors^ in which Mn acts as an isovalent 
magnetic impurity. The dependence on p in (Ga,Mn)As 
is expected to become very weak, however, when reach- 
ing the uncompensated limit or when further increasing 
hole density by non-magnetic acceptor co-doping. 



Microscopic TBA/CPA & mean-field approx. 
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FIG. 2: T c calculations within the microscopic TBA/CPA 
model: T c versus hole density (left panel), T c versus number 
of holes per MnGa (right panel) . 



D. Collective Mn-moments fluctuations 

The potential importance of correlated Mn-moment 
fluctuations on T c in (Ga,Mn)As can be illustrated by re- 
calling, within a simple parabolic band model, the RKKY 
(or Friedel) oscillations effect which occurs as a conse- 
quence of the 2k p anomaly in the wavevector dependent 
susceptibility of the hole systemiiM^ In this picture, the 
sign of the indirect kinetic- exchange MnQa-MnQa cou- 
pling fluctuates as cos(2/c^ d), where d is the distance 
between Mnc a moments, and its amplitude decays as d 3 . 
We can estimate the average Mnca-MnQa separation in a 
(Ga,Mn)As random alloy as d = 2(3/4nN Mn ) 1/3 - If the 



spin-orbit interaction and band-warping are neglected, 
the top of the valence band is formed by six degenerate 
parabolic bands with kp — {-K 2 p) 1 ^' i . For uncompen- 
sated (Ga,Mn)As systems (p = Nmu), we then obtain 
cos(2kpd) k — 1 which means that role of these fluctua- 
tions cannot be generally discarded. In realistic valence 
bands, as we see below, the fluctuations is suppressed 
due to non-parabolic and anisotropic dispersions of the 
heavy- and light-hole bands and due to the strong spin- 
orbit coupling. 

On a more quantitative level, we can establish the 
range of reliability and estimate corrections to the 
mean-field theory in (Ga,Mn)As by accounting for the 
suppression of the Curie temperature within quan- 
tum theory of long-wavelength spin-waves in the semi- 
phenomenological virtual-crystal model. We note that 
a qualitatively similar picture is obtained using Monte 
Carlo simulations which treat Mn-moments as classi- 
cal variables and account for positional disorderiiiMiiS 
Isotropic ferromagnets have spin-wave Goldstone collec- 
tive modes whose energies vanish at long wavelengths, 

n k = Dk 2 + 0(k 4 ) , (9) 

where k is the wavevector of the mode. Spin-orbit cou- 
pling breaks rotational symmetry and leads to a finite 
gap. According to numerical studies^ this gap is small 
however, much smaller than ksT c for example, and plays 
a role in magnetic fluctuations only at very low temper- 
atures. Spin-wave excitations reduce the total spin by 
one, at an energy cost that is, at least at long wave- 
lengths, much smaller than the mean-field value, Hmf- 
The importance of these correlated spin excitations, ne- 
glected by mean-field theory, can be judged by evaluating 
an approximate T c bound based on the following argu- 
ment which uses a Debye-like model for the magnetic 
excitation spectrum. When spin-wave interactions are 
neglected, the magnetization vanishes at the tempera- 
ture where the number of excited spin waves equals the 
total spin of the ground state: 

N Mn S = / ° dk k 2 n(fl k ) , (10) 
2tH Jo 

where n(Qk) is the Bose occupation number and the De- 
bye cutoff, kp, = (G^NMn) 1 ' 3 ■ It follows that the ferro- 
magnetic transition temperature cannot exceed 

29+1 

k B T c = —I—k 2 D D(T c ) . (11) 

In applying this formula to estimate T c we have approx- 
imated the temperature dependence of the spin stiffness 
by 

D(T) = D (S)(T)/S , (12) 

where Do is the zero-temperature stiffnesses and 
(S)(T) is the mean- field Mn polarization 2 - 5 at a tem- 
perature T. If the difference between T c obtained from 
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the self-consistent solution of Eqs. and l|12fl and the 
mean-field Curie temperature in Eq. Q is large, the typ- 
ical local valence-band carrier polarization will remain 
finite above the critical temperature and ferromagnetism 
will disappear only because of the loss of long-range spa- 
tial order, the usual circumstance for transition metal 
ferromagnetism for example. 

In discussing corrections to mean-field-theory T c esti- 
mates, we compare spin-stiffness results obtained with 
the simple two-band and realistic six-band models. De- 
tails on the formalism used to calculate Do can be found 
in Refs. I33l54l We find that the zero-temperature spin 
stiffness is always much larger in the six-band model. 
For (Ga,Mn)As, the two-band model underestimates D 
by a factor of ^10-30 over the range of hole densities 
considered. Furthermore, the trend is different: in the 
two-band model the stiffness decreases with increasing 
density, while for the six-band description the initial in- 
crease is followed by a saturation. Even in the limit of 
low carrier concentrations, it is not only the (heavy-hole) 
mass of the lowest band which is important for the spin 
stiffness. In the realistic band model, heavy-holes have 
their spin and orbital angular momenta aligned approxi- 
mately along the direction of the Bloch wavevector. Ex- 
change interactions with Mn spins mix the heavy holes 
with more dispersive light holes. The calculations show 
that heavy-light mixing is responsible for the relatively 
large spin stiffnesses. Crudely, the large mass heavy-hole 
band dominates the spin susceptibility and enables local 
magnetic order at high temperatures, while the dispersive 
light-hole band dominates the spin stiffness and enables 
long-range magnetic order. The multi-band character of 
the semiconductor valence band plays an important role 
in the ferromagnetism of these materials. 

Blue symbols in Fig. ^ summarize critical temperature 
estimates that include both the Stoner enhancement of T c 
and the suppression due to spin- wave fluctuations. The 
data were calculated using the six-band Kohn-Luttinger 
model for hole densities up to one hole per MnQ a and 
MnQa concentrations x — 5,8 and 10%. Given the qual- 
itative nature of these T c estimates we can conclude that 
T c will remain roughly proportional to x even at large 
dopings. The suppression of T c due to spin- waves in- 
creases with increasing hole density relative to the local 
moment concentration, resulting in saturation of the crit- 
ical temperature with increasing p at about 50% compen- 
sation. 



E. MriGa and Mni partial concentrations 

In the previous sections we have considered Mn to 
occupy only the Ga substitutional positions and found 
that T c should increase linearly with the concentration 
of MnQa local moments. In real (Ga,Mn)As materials a 
fraction of Mn is incorporated during the growth in in- 
terstitial positions. These donor impurities are likely to 
form pairs with MnQ a acceptors in as-grown systems with 



approximately zero net moment resulting in an ef- 

fective free local-moment doping x e ff = x s —Xi. Here x s 
and partial concentrations of substitutional and 

interstitial Mn, respectively. Although Mni can be re- 
moved by low-temperature annealing, a e // will remain 
smaller than the total nominal Mn doping. The MnQ a 
doping efficiency is, therefore, one of the key parame- 
ters that may limit maximum T c achieved in (Ga,Mn)As 
epilayers. 

In this Section, we calculate cohesion energy E c (x s , Xi) 
as a function of the partial concentrations x s and Xi 
and use it to determine the dependence of x s and Xi 
on the total Mn doping in as-grown materials. We define 
E c (x s ,Xi) as a difference of the crystal energy per unit 
cell and a properly weighted sum of energies of isolated 
constituent atoms. The cohesion energy is not very sen- 
sitive to the details of the electronic structure and can 
be calculated with a reasonable accuracy, for example by 
using the microscopic TBA model. Note that the growth 
kinetics calculations? 5 identified adsorption pathways for 
Mni formation in (Ga,Mn)As epilayers. Our equilibrium 
consideration provide, as seen in Section IIII Bl a very 
good estimate for the fraction of Mn impurities incorpo- 
rated in interstitial positions. 

The partial Mn concentrations x s and xi can be ob- 
tained by minimizing E c (x s ,Xi) at fixed Mn concentra- 
tion X - — - X £ I Xi with respect to either x s or Xi . Formally 
the condition for a dynamical equilibrium between the 
two positions of Mn has a form 

dE c (x s ,Xj) _ dE c (x s ,Xj) _ 
dx s dxi 

It was recently shown 5 ^ that the partial derivatives of 
the cohesion energy E c (x 8 , x$) with respect to x s and x$ 
represent formation energies F s and Fi of Mnc a and Mni 
impurities, respectively, assuming that the atomic reser- 
voir is formed by neutral isolated atoms. The equilibrium 
distribution of Mnc a and Mni is therefore reached when 

F s (x a ,Xi) = Fi(x a ,Xi), (14) 

as expected also from the growth point of view. Par- 
tial concentrations x St i of Mn can be obtained by solv- 
ing Eq. H13|) together with the condition < x St i < x, 

Xg \ X ^ " — - X, 

The left inset of Fig. [31 summarizes the compositional 
dependence of the cohesion energy in (Ga,Mn)As with 
both MnQa and Mni impurities. We consider several val- 
ues of x and plot E c (x s ,x — x s ) vs x s . Although the 
changes of the cohesion energy due to the incorporation 
of Mn are small, a systematic linear shift of the mini- 
mum of E c with increasing x is clearly visible. Corre- 
spondingly the partial concentration of x s ^ are expected 
to increase with increasing x. For x > 1.5% we obtain 
x s « 0.8a: and Xi « 0.2a;, in good agreement with the 
density-functional results^. 

The linear relations between x s , Xi, and x reflect the 
fact that the difference of the formation energies of MnQ a 
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F. Hole density and Hall coefficient 




FIG. 3: Main panel: Theoretical equilibrium partial con- 
centrations of substitutional MnGa (red line) and interstitial 
Mni (blue line) impurities. Right inset: Formation energies of 
MnGa and Mni as a function of total Mn concentration. Left 
inset: Cohesion energy as a function of substitutional MnGa 
concentration at several fixed total Mn concentrations. 



and Mni impurities (see right inset of Fig. EJ) can be, up 
to x = 10%, approximated by a linear function of x s and 

Xi-, 



&(x s ,Xi) = F s (x s ,Xi)-Fi(x s ,Xi) « -0.1+5. 9x s -15.1Xi(eV) 

(15) 

This relation allows us to interpret the theoretical dis- 
tribution of Mn atoms between substitutional and inter- 
stitial sites. For x < 1.5%, MnQ a has a lower formation 
energy than Mni and Mn atoms tend to occupy substi- 
tutional positions. At x w 1.5%, A.(x s ,Xi) approaches 
zero and both MnQ a and Mni are formed with a similar 
probability, as shown in Fig. |3J 

We note that both MnQ a and Mni positions remain 
metastable in the whole concentration range shown in 
Fig. and that our results correspond to the as-grown 
rather than to the annealed materials. During the 
growth, the formation energies (namely A.(x s ,Xi)) con- 
trol incorporation of Mn atoms assuming that the total 
amount of Mn in the material is related to a sufficiently 
high chemical potential in the Mn source. The annealing 
processes, on the other hand, do not depend on forma- 
tion energies but rather on energy barriers surrounding 
individual metastable positions of Mn in the lattice. The 
barriers are larger for MnQa^' 5 " so that the post-grown 
low temperature annealing can be used to remove Mni 
without changing the number of MnQ a significantly. 



As discussed above, the level of compensation is one 
of the key parameters that determines Curie temper- 
atures in (Ga,Mn)As. In this paper, as well as in a 
number of other experimental works, hole densities are 
obtained from Hall measurements. In order to esti- 
mate the uncertainty of this experimental technique we 
analyze in this section theoretical Hall factors, Th = 
(Pxy — Pxy,o)/(B/ep), in ferromagnetic (Ga,Mn)As epi- 
layers. Here p xy ,o is the Hall resistivity at field B = 
which can be non-zero due to the anomalous Hall effect. 

Detailed microscopic calculations in non-magnetic p- 
type GaAs with hole densities p ~ 10 17 — 10 20 cm' 3 
showed that r# can vary between 0.87 and 1.75, depend- 
ing on doping, scattering mechanisms, and on the level 
on which the complexity of the GaAs valence band is 
modelled^ 7 " Here we focus on estimating the effect on 
rs of the spin-splitting of the valence band and of the 
anomalous Hall term that is particularly large in ferro- 
magnetic (Ga,Mn)As. 

The calculations are based on numerical evaluation 
of the Kubo formula at finite magnetic fields. We as- 
sume band- and wavevector-independent quasiparticle 
lifetimes for simplicity. It is essential for our analy- 
sis to allow for both intra-band and inter-band tran- 
sitions. At zero magnetic field, the inter-band transi- 
tions between SO-coupled, spin-split bands give rise to 
the anomalous Hall effect (AHE), i.e., to a non-zero p xy 
that is proportional to the magnetization^ On the other 
hand, the ordinary Hall resistance which is proportional 
to B arises, within the simple single-band model, from 
intra-band transitions between adjacent Landau levels. 
The Kubo formula that includes both intra-band and 
inter-band transitions allows us to capture simultane- 
ously the anomalous and ordinary Hall effects in the com- 
plex (Ga,Mn)As valence bands. 

Many of the qualitative aspects of the numerical data 
shown in Figs. 01 and [S] can be explained using a sim- 
ple model of a conductor with two parabolic uncou- 
pled bands. Note that the typical scattering rate in 
(Ga,Mn)As epilayers is H/t ~ 100 meV and the cyclotron 
energy at B = 5 T is fku ^ 1 meV, i.e., the system is in 
the strong scattering limit, ujt « 1. In this limit, the 
two band model gives resistivities: 



1 



1 



Pxy 



&xx,l + <?xx,2 fO,l+Co,2 
°xy,\ + &xy,2 



{Pxx,\ + Oxx.IlY 



B ! + >B 
epi (1 + 21 ml)2 ~ ep 



(16) 



where the indices 1 and 2 correspond to the 1st and 2nd 
band respectively, the total density p = P\ + P2, and the 
zero-field conductivity 00 = e 2 Tp/m*. Eq. Ijl6(l suggests 
that in the strong scattering limit the multi-band nature 
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of the hole states in (Ga,Mn)As should not result in a 
strong longitudinal magnetoresistance. This observation 
is consistent with the measured weak dependence of p xx 
on B for magnetic fields at which magnetization in the 
(Ga,Mn)As ferromagnet is saturated^ 

The simple two-band model also suggests that the Hall 
factor, rn, is larger than one in multi-band systems with 
different dispersions of individual bands. Indeed, for un- 
coupled valence bands, i.e. when accounting for intra- 
band transitions only, the numerical Hall factors in the 
top panels of Figs. 0] and are larger than 1 and inde- 
pendent of r as also suggested by Eq. (|Td|| . The suppres- 
sion of rn when SO-coupling is turned on, shown in the 
same graphs, results partly from depopulation of the an- 
gular momentum j = 1/2 split-off bands. In addition to 
this "two-band model" like effect, the inter-Landau- level 
matrix elements are reduced due to SO-coupling since 
the spinor part of the eigenfunctions now varies with the 
Landau level index. In ferromagnetic Gai-^Mn^As the 
bands are spin-split and higher bands depopulated as x 
increases. In terms of ru, this effect competes with the 
increase of the inter-Landau-level matrix elements since 
the spinors are now more closely aligned within a band 
due to the exchange field produced by the polarized Mn 
moments. Increasing x can therefore lead to both de- 
crease or increase of rn depending on other parameters 
of the (Ga,Mn)As, such as the hole density (compare top 
right panels of Figs. 0] and [5]). 

The inter-band transitions result in a more single-band 
like character of the system, i.e. rn is reduced, and the 
slope of the p xy {B) curve now depends more strongly 
on r. Although the anomalous and ordinary Hall ef- 
fect contributions to p xy cannot be simply decoupled, the 
comparison of numerical data in the four panels and the 
inset in Fig. [5] confirms the usual assumption that the 
anomalous Hall effect produces a field-independent off- 
set proportional to magnetization and p xx . The compar- 
ison also suggests that after subtracting p xy (B = 0), rn 
can be used to determine the hole density in (Ga,Mn)As 
with accuracy that is better than in non-magnetic GaAs 
with comparable hole densities. For typical hole and Mn 
densities in experimental (Ga,Mn)As epilayers we esti- 
mate the error of the Hall measurement of p to be within 
±20%. 



III. EXPERIMENT 



Measured Curie temperatures and hole 
densities 



A series of (Ga,Mn)As thin films with varying Mn con- 
tent were grown by LT-MBE epitaxy using As2. The 
layer structure is 25 or 50nm (Ga,Mn)As / 50nm low 
temperature GaAs / lOOnm high temperature (580°C) 
GaAs / GaAs(lOO) substrate. For a given Mn content, 
the growth temperature of the (Ga,Mn)As layer and the 
GaAs buffer is chosen in order to minimize As antisite 
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FIG. 4: Theoretical Hall factors for p = 0.8 nm" 3 ; h/r = 
50 meV (dashed lines), H/t = 150 meV (solid lines). Top 
panels: only intra-band transitions are taken into account. 
Bottom panels: intra- and inter-band transitions are taken 
into account. Left panels: GaAs (a; = 0); zero SO-coupling 
(red lines), Aso = 341 meV (blue lines). Right panels: 
(Ga,Mn)As with MnGa concentration 4% (green lines), 8% 
(brown lines). p xy — in all panels except for the bottom 
left panel where p xy (B = 0) ^ due to the anomalous Hall 
effect. 



densities while maintaining two-dimensional growth and 
preventing phase segregation. We find that the growth 
temperature must be decreased as the Mn concentration 
is increased: for the lowest Mn content the growth tem- 
perature was ~ 300°C, for the highest it was ~ 180°C. 
Full details of the growth are presented elsewhere ^SiSi 

The Mn content was controlled by varying the Mn/Ga 
incident flux ratio, measured in — situ and calibrated us- 
ing SIMS measurements on lam thick (Ga,Mn)As films, 
grown under otherwise identical conditions to the sam- 
ples considered here. A detailed comparison of the results 
of a number of different calibration techniques, presented 
in Ref. |62| allows us to assign an uncertainty of ±10% to 
the quoted Mn concentrations. However, it should be 
noted that the SIMS measurements yield the total vol- 
ume density of Mn in the (Ga,Mn)As films, and not the 
fraction of Ga substituted by Mn. This is important as it 
is expected that a fraction of the Mn will be incorporated 
on interstitial as well as substitutional sites^i. We define 
the Mn concentration, the total Mn volume density 
relative to the volume density of Ga in GaAs. 

Hall bar structures, of width 200 pm and length 1 
mm, were fabricated from the (Ga,Mn)As samples using 
photolithography. Simultaneous magnetoresistance and 
Hall effect measurements were performed using standard 
low-frequency ac techniques, in order to extract both the 
Curie temperature T c and the hole density p, as detailed 
below. Magnetic fields of up to ±0.7 T and ±16.5 T 
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FIG. 5: Theoretical Hall factors for p = 0.2 nm -3 ; same line 
coding as in Fig. [I] Inset: Theoretical Hall curves showing 
the anomalous Hall effect contribution at B = 0. 



were used to obtain T c and p, respectively. Following 
these measurements, the samples were annealed in air 
at 190°C. The electrical resistance was monitored during 
annealing, and the anneal was halted when this appeared 
to have reached a minimum (after typically 50 to 150 
hours). The T c and p were then re- measured. 

Below T c , the Hall resistance R xy in (Ga,Mn)As is 
dominated by the anomalous Hall effect, with R xy ~ 
RaM z , where M z is the perpendicular component of the 
magnetization, and the coefficient Ra is roughly propor- 
tional to the square of the resistivity, p xx . Therefore, R xy 
I p xx gives a direct measurement of M z , which can be 
used to extract T c using Arrot plotsS The value of T c 
obtained depends only weakly on the precise dependence 
of Ra on p xx assumed, since p xx varies only slowly close 
to T c , while R xy varies rapidly. We are therefore able to 
obtain T c within an accuracy of ±1 K by this method^ 

T c obtained for the (Ga,Mn)As Hall bar samples be- 
fore and after annealing are shown versus x in Fig. [5] 
It can be seen that the low-temperature annealing pro- 
cedure results in a marked increase in T c as has been 
found previously^ Increases of T c by more than a fac- 
tor of two are possible. This effect becomes larger as 
the Mn concentration increases. Since the T c enhance- 
ment is associated with out-diffusion and passivation of 
interstitial Mnji this indicates that as the incident Mn 
flux is increased, an increasing fraction is incorporated 
on interstitial sites, as predicted in Section III El 

To obtain hole densities from R xy , it is necessary to 
separate the small normal Hall term from the much larger 
anomalous Hall term. Measurements were performed at 
0.3 K and in magnetic fields above 10 T, i.e., under con- 
ditions where the normal Hall term gives the dominant 
field-dependent contribution to R xy . Then, the measured 



FIG. 6: Experimental Curie temperature versus total Mn 
doping. T c is measured from the anomalous Hall effect, Mn 
doping by SIMS. Open symbols correspond to as-grown sam- 
ples, half-open symbols to as-grown samples with large charge 
compensation, and filled symbols to annealed samples. For 
clarity, error bars are shown only for the x=6.7% sample. 

R xy was fitted to (ap 2 + ruB), where p xx and B are the 
measured resistivity and magnetic field, and a and th 
are fit parameters. Finally, the hole density is obtained 
from m = l/(pew), where w is the (Ga,Mn)As layer 
thickness. From the detailed calculations described in 
Section III Fl we can ascribe an uncertainty of ±20% to 
the values of p obtained using this method. The mea- 
sured p for the (Ga,Mn)As Hall bar samples before and 
after annealing are shown versus x in Fig. d We see that 
annealing greatly increases p for large x. Data in the in- 
set of Fig. [7] discussed in detail below, show that within 
error the samples are uncompensated after annealing. 



B. Substitutional and interstitial Mn 

From the measured hole density p before and after an- 
nealing, and the total Mn density x, values can be ob- 
tained for the density of incorporated Mn occupying ac- 
ceptor substitutional and double donor interstitial lattice 
sites, x s and Xi. These are obtained using the following 
assumptions: i) the only contribution to the total Mn 
density determined by SIMS are from substitutional and 
interstitial Mn, i.e. x = x s + Xjj ii) the only source of 
compensation in the (Ga,Mn)As films are the interstitial 
Mn, which are double donors i.e. p = A/af c (x s — 2xi); 
iii) the low temperature annealing procedure affects only 
Xi, and not x s . The values of x s and Xi in the unan- 
nealed films obtained under these assumptions are shown 
in Fig. [SJ We find a remarkably good agreement between 
experiment and the theory data in Fig.|21and in ReflSSl 
As a consistency check, we show in the inset of Fig. [7| 
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FIG. 7: Experimental hole density versus total Mn doping. 
Hole density is measured by ordinary Hall effect. Same sym- 
bol coding is used as in Fig|S| 




FIG. 8: Experimental partial concentrations of MnGa (trian- 
gles) and Mni (circles) in as grown samples. Data show no 
saturation of MnGa with increasing total Mn doping consis- 
tent with theory expectation. 



the ratio of hole density to substitutional Mnc a density 
after annealing, as obtained under the above assump- 
tions. Within the experimental error we obtain one hole 
per substitutional MnQ a after annealing, that is, there is 
no significant compensation in the annealed (Ga,Mn)As 
films. This justifies our neglect of additional compensat- 
ing defects such as Asc a in determining x s and Xi. 



C. T c versus MnGa, effective Mnc a , and hole 
densities 

Since we obtain reasonably accurate values for T c , hole 
densities and the partial Mnc a and Mni concentrations 
for the set of samples considered here, we now attempt 
to assess on the basis of the experimental data the key 
factors determining T c and to compare the experimental 
results with the broad predictions of theory. 

In Fig. El T c was plotted against total Mn concentra- 
tion. Before annealing the T c values of samples with 
high compensation (samples with large compensation are 
indicated as half filled symbols in this and subsequent 
figures) do not increase significantly with increasing to- 
tal Mn density but a steady increase is recovered after 
annealing. In Fig. [5] T c is plotted against the substitu- 
tional MnGa concentration. The form of Figs. HO and [5] 
are broadly similar despite the different x-axes. We ex- 
pect, however, and will assume in the following discussion 
that any Mni donor present is attracted to a Mnc a ac- 
ceptor and that the pair couples antiferromagnetically. 37 
Then the effective uncompensated moment density will 
be Plotting T c against x e f / in Fig. EH 

reveals that for all the low compensation samples T c in- 
creases approximately linearly with x e f / but that as com- 
pensation, (1 — paf c /4x e ff), increases above ~ 40% the 
measured T c values fall increasingly far below this linear 
trend. 

If T c is plotted against hole density, as is done in the 
inset of Fig. ^] it is found to increase monotonically. 
However, this is primarily due to the increase in hole 
density with x e f / . The main plot in Fig. 1111 shows that 
T c /x e f f is almost independent of hole density except for 
the case of the high compensation samples which clearly 
stand out as showing different behavior. To compare 
with the predictions of Section^we finally plot T c /x e ff 
against p/x e ff in Fig. ^| All experimental points in 
this plot show a common T c trend and the magnitudes of 
the experimental and calculated T c /x e ft are comparable. 
Further confirmation of the theoretical picture is seen 
from the very weak experimental dependence of T c /x e ff 
on p/x e ff for low compensation and the relatively rapid 
fall of T c /x e ff with decreasing p/x e ff for compensations 
of ~ 40% or larger. 

As a consistency check for considering x e f / as the den- 
sity of local MnQa moments participating in the ordered 
ferromagnetic state, magnetization data are shown in in- 
sets of Figs. 191 and ITU1 Magnetizations were determined 
by superconducting quantum interference device magne- 
tometry, at a sample temperature of 5K, and using an 
external field of 0.3 T to overcome in-plane anisotropy 
fields. The charge and moment compensation after an- 
nealing is not significant for our samples and the moment 
per x s or x e ff is within error 4-4.5 [Ib- This corresponds 
well to the 5 fiB contribution of the S=5/2 local MnQ a 
moment and an approximately (-0.5)-(-0.8) [Ib contribu- 
tion of the antiferromagnetically coupled valence band 
hole 65 in collinear (Ga,Mn)As ferromagnets. In the inset 
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FIG. 9: Experimental T c versus MriGa concentration, x s (see 
text for definition of x s ). Magnetization per x s is shown in 
the inset. 



FIG. 11: Experimental T c /x e ff versus hole density. T c /x e // 
is nearly independent of hole density except in highly com- 
pensated samples. Inset: T c versus hole density. 
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FIG. 10: Experimental T c versus effective MnGa concentra- 
tion, x e ff (see text for definition of x e ff). Magnetization per 
Xeff is shown in the inset. 



of Fig. [5] we see that the measured moment per x s are 
all below 4 /is for the compensated samples. Including 
the effects of the Mni-Mnca antiferromagnetic coupling 
by considering the moment per x e f / reveals again values 
around 4.5 [Ib- Our conclusion therefore is that if we 
assume no significant frustration in our samples and ac- 
count for the antiferromagnetic Mni-MnQ a coupling, our 
extensive set of T c , hole density, Mn density, and mag- 
netization data brings up a clear common picture of T c 
trends in the 14 different (Ga,Mn)As ferromagnetic semi- 
conductors we have studied, that is consistent with the 
theory predictions summarized in Section ITU 



FIG. 12: Experimental T c /x e ff versus hole density relative 
to effective concentration of MnGa moments. Deviation from 
linear dependence on x e ff are seen only for high compensa- 
tions (1 — paf c /4:Xeff > 40%) in agreement with theory. For 
weakly compensated samples T c shows no signs of saturation 
with increasing x e fj. 



IV. DISCUSSION 

The preceding considerations of the factors determin- 
ing T c in (Ga,Mn)As lead us to conclude that there are 
no fundamental physics barriers to achieving room tem- 
perature ferromagnctism in this system. Experimental 
results for T c in samples in which compensating defects 
other than interstitial Mn have been reduced to very low 
levels have been shown to be in good agreement with 
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theoretical expectations. Moment compensation by in- 
terstitial Mnj impurities becomes increasingly important 
as the concentration of total Mn is increased. However, 
for the range of total Mn concentrations considered ex- 
perimentally we find that the level of substitutional Mnc a 
continues to increase with x. Furthermore low tempera- 
ture post growth annealing is found to effectively remove 
Mni in thin film samples even at large x, leading to ma- 
terial which within experimental error is both charge and 
moment uncompensated. Most importantly for samples 
in which the charge compensation is less than ~ 40% we 
find theoretically and experimentally that T c increases 
approximately linearly with effective concentration, x e ff, 
of MnQa whose moments are not compensated by near- 
neighbor Mni impurities. We have not observed any signs 
of saturation in this trend in the studied (Ga,Mn)As di- 
luted magnetic semiconductors. It should be noted that 
our maximum x e ff is only 4.2% in the as grown sample 
and 6.2% after annealing for a total Mn concentration 
x = 9%. Hence the modest T c 's observed so far. Achiev- 
ing T c values close to room temperature in (Ga,Mn)As, 
which we expect to occur for x e ff w 10% is essentially 
a technological issue, albeit a very challenging one. In 
the remaining paragraphs of this section we discuss these 
challenges in more detail. 

Low temperature MBE growth is used to achieve levels 
of Mn incorporation in (Ga,Mn)As far in excess of the 
equilibrium solubility level. When growing (Ga,Mn)As 
with Mn concentrations of several percent it is known 
that the Mn tends to accumulate on the surfac o 60 i 61 in a 
similar way to all high vapor pressure dopants in GaAs 
and to the higher vapor pressure species, e.g. In in In- 
GaAs. For homogenous Mn incorporation during con- 
tinuous growth, a surface Mn concentration is required 
that is higher than the bulk concentration. For a given 
Mn concentration this density gradient is temperature 
dependent, increasing with increasing temperature. This 
leads to an upper temperature limit for successful growth 
when the Mn surface concentration approaches a signifi- 
cant proportion of a monolayer, after which point surface 
clustering of Mn occurs, frustrating the growthiSfliSi Fur- 
thermore, higher Mn fluxes require lower growth temper- 
atures. 

The pursuit of higher Curie temperatures has driven 
growth efforts to very low temperatures compared with 
conventional MBE of GaAs. In this regime (^200 - 
250° C) significant levels of compensating defects such as 
Asca and vacancies usually occur in GaAsj£& The density 
of Asca defects can be reduced by close to stoichiomet- 
ric growth with As2 requiring very precise control over 
the As flux. 

Apart from precise control over the stoichiometry and 
the attendant requirement for flux stability, a major tech- 
nical difficulty arises from the measurement and control 
of the growth temperature. In order to measure sub- 
strate temperature, most MBE machines in use today 
employ a thermocouple heated by radiation from the sub- 
strate or substrate holder. At normal growth tempera- 



tures (~ 580° C) the radiant flux from the substrate is 
high and the relationship between substrate and thermo- 
couple is repeatable with a short time constant, allow- 
ing for good temperature stability and control. At low 
temperatures, however, the radiant flux between the sub- 
strate and thermocouple is low, leading to a heightened 
sensitivity to local conditions such as holder emissivity, 
radiant heat from the metal sources, shutter transients 
etc., and also long time constants. This significantly in- 
creases the error in the temperature measurement as well 
as the likelihood of temperature spikes and drift as shut- 
ters are opened and growth proceeds. In MBE, optical 
pyrometers are ubiquitous as secondary temperature cali- 
bration devices but most cannot read accurately at these 
low temperatures and in many common configurations 
suffer from potential inaccuracies due to reflection off the 
Knudsen cells if used during growth. 

It is desirable to grow at as high temperature as possi- 
ble for a given Mn flux, while maintaining 2D growth and 
avoiding Mn clustering. However, with such large errors 
and potential temperature drift, growers tend to err to- 
wards lower than ideal temperatures in order to sustain 
the growth. To explore fully the parameter space, ef- 
fort should be directed towards improving the control of 
both metal fluxes and substrate temperature. This will 
maximize the chances of increasing the doping towards 
the 10% MnQa, required for room temperature ferromag- 
netism. The increases in T c achieved in the last few years 
lead us to believe that higher transition temperatures will 
be obtained using conventional MBE. Growth interrupt 
strategies such as migration enhanced epitaxy (MEE) 68 
may have advantages over conventional MBE for the in- 
corporation of higher levels of substitutional Mn however 
they will be especially sensitive to poor temperature sta- 
bility and shutter transients and so will require even more 
precise temperature control. 



V. CONCLUSION 

Based on the broad agreement between theoretical and 
experimental Curie temperature trends in (Ga,Mn)As 
with Mn concentrations larger than 1.5% we can outline 
the following strategies for achieving room temperature 
ferromagnctism in this semiconductor: 

(i) T c increases linearly with the concentration, x e ff, 
of local MnQa moments participating in the ordered fer- 
romagnetic state. Room temperature ferromagnetism 
should be achieved at x e ff w 10%. Interstitial Mni im- 
purities reduce the number of these ordered MnQa mo- 
ments due to the strong antiferromagnetic Mnca-Mni 
near-neighbor coupling. Mni, however, can be efficiently 
removed by post-growth annealing, (ii) Equilibrium con- 
siderations, confirmed experimentally in samples with 
MnQa concentrations up to 6.2%, suggest that there is 
no fundamental physics barrier for increasing MnQa con- 
centration to and beyond 10%. A very precise control 
over the growth temperature and stoichiometry is, how- 
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ever, required for maintaing the 2D growth mode of the 
high quality, uniform (Ga,Mn)As materials, (iii) Ferro- 
magnetic coupling between the ordered local MnQ a mo- 
ments is mediated by itinerant holes. For charge com- 
pensations (1 — paf c /4:X e ff) > 40%, the Curie temper- 
ature falls down with decreasing p. At compensations 
smaller than ~ 40%, however, T c is almost independent 
of the hole density. A modest charge compensation is, 
therefore, not an important limiting factor in the search 
of high Curie temperature ferromagnetic semiconductors 
based on (Ga,Mn)As and may be desirable to maximize 
the possibilities for doping and gate control of ferromag- 
nctism. 
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